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Abstract 

The possibility of detecting Gamma-Ray Bursts (GRBs) in TeV energy 
range using large area muon detectors like AMANDA and Lake Baikal detector 
is examined. These detectors can detect TeV energy photons by detecting the 
secondary muons created by the TeV photons in the Earth's atmosphere. We 
calculate the expected number of muons and the signal to square root of noise 
ratios in these detectors due to TeV gamma-rays from individual GRBs for 
various assumption on their luminosity, distance from the observer (redshift), 
gamma-ray integral spectral index, maximum energy cutoff of the photon spec- 
trum, and duration, including the effect of the absorption of TeV photons in the 
intergalactic infrared radiation background. We also calculate the expected rate 
of detectable TeV-GRB events in these detectors using a recent determination of 
the luminosity and redshift distributions of the GRBs in the Universe. For rea- 
sonable ranges of values of various parameters, we find about 1 event in 20 years 
in AMANDA (and a similar number for BAIKAL), while the event rate can be 
significantly larger (by factors of 10 or more depending on the area of the detec- 
tor) in the proposed next generation detectors such as ICECUBE. Although, 
for the specific forms of the luminosity- and redshift distributions assumed, the 
average rate of expected detection is low, occasional nearby {z < 0.1), high lu- 
minosity ( > 10^^ erg/ sec), long duration ( > lO's of seconds), and sufficiently 
hard spectrum (gamma-ray integral spectral index < 1) GRBs can, however, be 
detected even by AMANDA. Detection (or even non-detection) of TeV photons 
from GRBs in coincidence with satellite-borne detectors (which detect mainly 
sub-MeV photons) will be able to provide important new insights into the char- 
acteristics of GRBs and their emission mechanisms and, in addition, provide 
limits on the strength and spectrum of the intergalactic infrared background 
which affects the propagation of TeV photons from cosmological sources. 
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1 Introduction 



Gamma Ray Bursts (GRBs) are among the most powerful astrophysical phenomena 
in the Universe; see e.g., Ref. |I| for a recent review. Because of the hmited sizes 
of the satelhte-borne GRB detectors (BATSE, BeppoSAX, for example), GRBs have 
been detected mostly in the sub-MeV energy region where the photon number flux 
is sufficiently large for typical power-law photon spectra of GRBs. However, GRB 
emission extending to ~10 GeV is known as in the case of the famous long-duration 
burst GRB 940217 detected by the EGRET instrument on board the Compton 
Gamma Ray Observatory (CGRO), suggesting the possibility that GRBs may also 
emit much higher energy photons, perhaps even extending to TeV energies as in some 
highly energetic Active Galactic Nuclei (AGN) of the "blazar" class. For power-law 
spectra falling with energy, the photon number flux at TeV energies may be too low 
for TeV photons from GRBs to be detected by the satellite-borne detectors. However, 
ground-based detectors can in principle detect TeV photons from GRBs by detecting 
the secondary particles (electrons, muons) of the atmospheric showers generated by 
these photons in the Earth's atmosphere. A variety of techniques (see, e.g., for 
a review) allow determination of the energy and direction of the shower-initiating 
primary photon with good accuracy. 

Indeed, three major ground-based gamma ray detectors, the Tibet air shower array [Q, 
the HEGRA-AIROBICC Cherenkov array |0 and the Milagro water-Cherenkov de- 
tector have independently claimed evidence for possible TeV 7-ray emission from 
sources in directional and temporal coincidence with some GRBs detected by BATSE. 
The estimated energy in TeV photons in each case has been found to be about 1 to 
2 orders of magnitude larger than the corresponding sub-MeV energies measured by 
BATSE. Since TeV photons are efficiently absorbed in the intergalactic infrared (IR) 
background due to pair production 0, only relatively close by (i.e., low redshift) 
GRBs, for which the absorption due to IR background is insignificant, can be ob- 
served at TeV energies, which would explain the fact that only a few, not all, of the 
BATSE-detected GRBs in the fields of view of the individual ground detectors have 
been claimed to be detected at TeV energies. 

While more GRBs would have to be detected in TeV energies before any firm con- 
clusions can be drawn, the above discussions already point to the possibility that not 
only may the spectra of some, if not all, GRBs extend to TeV energies, but that the 
energy emitted in TeV photons may even constitute the dominant part of the total 
energy emitted by individual GRBs. This would happen, for example, if the pho- 
ton spectrum, extending to TeV region, is hard with a differential power-law index 
7 < 2. In such situations, while the total number of photons in the entire spectrum 
is dominated by those at the lower cutoff energy of the spectrum, the total energy is 
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dominated by the upper cutoff energy. In other words, it is possible that the total 
burst energies of individual GRBs calculated from the measured fiuences in the keV 
- MeV energies by the satellite detectors may only be a small fraction of the actual 
total energies of the individual bursts, most of which could be carried away by the 
relatively few TeV photons rather than by the abundant sub-MeV photons. 

Indeed, photon differential spectrum with power-law index 7 < 2 would be expected 
in most situations if the radiation arises from synchrotron radiation of power-law 
distributed particles (electrons or protons): Recall that the power-law index of the 
synchrotron photon differential spectrum is (p+ l)/2 where p is the power-law index 
of differential number spectrum of the particles (electrons or protons). Thus, one can 
expect 7 < 2 as long as p < 3, which is typically the case for particles accelerated in 
shocks. 

On the theoretical side, it has been suggested P, ^ that if protons are accelerated to 
energies above ~ 10^° eV in GRBs as envisaged in the scenario |Ty] in which GRBs are 
the sources of the ultrahigh energy cosmic rays extending to energies above 10^° eV, 
then the synchrotron radiation of those ultrahigh energy protons in the magnetic field 
within the GRB would produce TeV photons. Thus, while the synchrotron radiation 
of electrons produce the keV-MeV flux, the TeV flux would be produced by syn- 
chrotron radiation of the protons with energy above 10^° eV. Whether or not a GRB 
would be bright in TeV would, in this scenario, be determined by the efficiency of 
energy transfer from protons to electrons within the GRB [^]. It has also been sug- 
gested [0 |l^ that while the TeV photons emitted by GRBs at large redshifts would 
be absorbed through e+e~ pair production on the intergalactic infrared background, 
the resulting electromagnetic cascades initiated by the produced pairs could produce 
the observed extragalactic diffuse gamma ray background in the GeV energy region. 
It is thus clear that confirmation of TeV gamma ray emission from GRBs would have 
major implications not only for the physics and astrophysics of GRBs, but also for 
the strength and redshift evolution of the intergalactic infra-red background. 

High energy primary gamma rays above a few hundred GeV create "air-showers" of 
secondary particles by interacting with Earth's atmosphere. Among the secondary 
particles are muons which can be detected, and the energy and direction of the shower 
initiating parent photon reconstructed (to a directional accuracy of typically a degree), 
by using the relatively shallow underground muon detectors such as AMANDA [ [I3|| , 
Lake Baikal detector [|14| and the proposed ICECUBE [|1^] , and surface muon detector 



such as MIL AGRO . The first three are facilities primarily for detecting high energy 
neutrinos from cosmic sources — the neutrinos are detected by detecting the muons 
produced by the neutrinos in ice (AMANDA and ICECUBE) or water (Lake Baikal) 
— while MILAGRO is a water-Cherenkov detector of "all" particles (including muons) 
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primarily for detecting gamma rays above a few hundred GeV from astronomical 
sources. 



The advantages of using these muon detectors over the conventional air-Cherenkov 



telescopes for doing TeV gamma ray astronomy have been expounded in Refs. [|T6|, [17 
Compared to air-Cherenkov telescopes, the muon detectors cover much larger fraction 
of the sky with large duty cycles. For example, the AMANDA detector at the South 
Pole covers more than a quarter of the sky with essentially 100% efficiency. Being at 
relatively shallow depths — of order 1 km — the "neutrino" detectors are sensitive to 
muons with energies of a few hundred GeV and thereby to parent photons of energy 
of order a TeV and above. MILAGRO detector's muon detection threshold is smaller 
— about 1.5 GeV — because of its location on the surface, which makes it sensitive 
to even lower energy gamma rays. 



As pointed out in Ref. IjT^], these muon detectors are particularly suited for detecting 
TeV gamma rays from transient sources like GRBs — the otherwise large background 
of down-going cosmic-ray induced atmospheric muons can be significantly reduced 
because, with the information on the time and duration of a burst provided by satellite 
observation, the background is integrated only over the relatively short duration of a 
typical GRB which is of order few seconds. 



In this paper we calculate, following Ref. fl^, the expected number of muons and the 



signal to square root of noise ratios in AMANDA and Lake Baikal detectors due to 
TeV gamma-rays from individual GRBs for various assumption on their luminosity, 
distance from the observer (redshift), gamma-ray integral spectral index, maximum 
energy cutoff of the photon spectrum, and duration, including the effect of the absorp- 
tion of TeV photons in the intergalactic infrared radiation background. We carefully 
take into account all relevant cosmological effects (in a cosmological constant domi- 
nated present Universe, see below for details). Our results for number of muons for 
individual GRBs differ significantly (especially in the case when absorption in the 



intergalactic IR background is taken into account) from those obtained in Ref. ||T7| 
for the same set of GRB parameters. We also calculate the expected rate of de- 
tectable TeV-GRB events in AMANDA-like detectors using a recent determination 
of the luminosity and redshift distributions of the GRBs in the Universe. 

This paper is organized as follows: In Sec. 2 we discuss the parametrization of the 
TeV emission from the GRBs and the formalism we have used to calculate the number 
of secondary muons produced in the atmosphere from these TeV gamma rays of the 
GRBs. In Sec. 3 the results of our calculation of the number of muons expected and 
the signal to square root of noise ratios due to TeV photons from individual GRBs as 
functions of various GRB parameters are discussed. The expected rate of detection of 
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TeV GRBs is calculated in Sec. 4 by using a recently determined luminosity function 
and redshift distribution of the GRBs in the Universe, with summary and conclusions 
presented in Sec. 5. 



2 Parametrizations of the Photon Spectrum of 
GRBs and the Secondary Muons in Gamma- Ray 
Induced Showers in Earth's Atmosphere 

2.1 Photon Spectrum of GRBs 

We shall assume the emission spectrum of the GRB to be constant during the emission 
time 6te, and represented by dN^/dE^ = KE~"~^ , giving the total number of photons 
emitted during the entire burst per unit energy at energy between a minimum and 
a maximum energy, -E-ymin and -E^max, respectively. Here a > is the integral spectral 
index and K isa. normalization constant, which are related to the total energy emitted 
from the burst, -Egrb, by 



from which we get 



Egrb = K E-'^dE,, (i; 



K = EoKB X <( + 1) (^-n-i' -^mii^T' " ^ 1 (2) 

[In (^^max/^7min)] ^ for « = 1 . 



Following Ref . [T^ , we shall use -E^min = 1 MeV since negligible amount of energy is 
emitted below that energy, and assume that the spectrum continues to TeV energies 
with the same integral spectral index a. Further, as discussed in the Introduction, 
we shall assume the photon spectrum to be sufficiently hard with differential spectral 
index 7 < 2, i.e., integral spectral index a < 1. 

The luminosity, L, and the total energy emitted, -Egrb, are related as 

Egrb = Lx6te. (3) 

For a GRB source located at a cosmological distance corresponding to a redshift z, 
the observed burst duration is St = {1 + z)6te. For our purpose, a GRB is completely 
specified by the parameters L, z, 6te, a and -E^max- 

The observed spectrum is determined by two factors: (a) the photons energies are 
redshifted by a factor (1 + z) due to the expansion of the universe, and (b) GeV 
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- TeV photons are absorbed in the intergalactic space due to e+e~ pair production 
on the intergalactic starhght and infrared background photons. The probabihty for a 
photon emitted with energy by a source at redshift z to reach the earth is e~'^'^^'"^\ 
where r is the optical depth for absorption due to the pair production process. We 
have used the optical depths given in [0 , where a parametrization has been given for 
r(i?^, z) for small z which, as we shall see, is appropriate for our purpose. 

The observed photon energy E^q and its energy at the source, E^, are related as 
E^Q = E^/{1 + z). Thus, at the top of the atmosphere, the observed cutoff energy of 
the source is -E^maxo = -E^max/ll + z). Assuming isotropic emission from the source, 
the total number of photons from a single GRB at redshift z per unit energy at energy 
E^Q that strike per unit area of the Earth at the top of the atmosphere is 

"^^^ ^ E-"-i(l + z)e-^(^-^), (4) 



dE^Q A7rr^{z) ^ 

where r{z) is the comoving radial coordinate distance of the source. For a spatially 
flat Universe (which we shall assume to be the case) with + Qm = 1, where Qa 
and Qm are, respectively, the contribution of the cosmological constant and matter to 
the energy density of the Universe in units of the critical energy density, 3Hq/8ttG, 
the radial coordinate distance r{z) is given by ^ 



Here c is the speed of light, Hq is the Hubble constant in the present epoch. In our 
calculations, we shall use = 0.7, VLrn = 0.3, and Hq = 65 km sec~^ Mpc~^. 

2.2 Number of Muons in a Photon-Induced Air Shower 

Photons of sufficiently high energy striking the earth's atmosphere interact with the 
air nuclei to produce pions; these photo-pions are the major source of muon production 
by photons in atmospheric air-showers, among other possible physical processes like 
direct pair production of muons by photons or photo-produced charm decays. The 
number of muons with energy above E^ in an air-shower produced by a single photon 
of energy E^q striking the atmosphere can be parametrized [1^, [1^, for E^j^ in the 
range 0.1 TeV to 1 TeV, by the formula 

^ X 2.14 X 10-^ 1 E^o 

N^iE,o, > E^) ~ /° , 6 

cost/ [E^/ COS 0) {Efj_/ COS 0) 

where 6 is the zenith angle of the photon source, and all energies are in TeV units. 
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The above parametrization is valid for E^q/E^ > 10 |T^, which is adequate for the 
purpose of calculating the muon rates in the AMANDA and Baikal detectors both of 
which have muon thresholds of the order of a few hundred GeV UTZ . 



The total number of muons with energy in excess of E^ in a detector of effective area 
A due to a single GRB at redshift z can then be written as 



A 



xO{2) 



dE. 



■70 



'7th 



dE. 



70 



{z)N^{E^o. > E, 



(7) 



where E^ytt ~ 10 x E^/ cos6 is the minimum photon energy needed to produce muons 
of energy E^^ in the atmosphere 6 being the zenith angle of the photon source. 



Using Eq. in (J^) we finally have 



A 



K{l + z)- 
A.Txr'^^z) 



xO(^) 



'7th 



e-(i^-,oi^+^MN^{E^,, > E,) 



where A' is a function of the GRB parameters L, z, 5te, a and iJ-ymax, and is given 
by Eqns. ^ and ^. 



3 TeV Photon Signal of GRBs in AMANDA and 
Lake Baikal Muon Detectors 

The threshold muon energy for detection of vertical muons in the AMANDA detec- 
tor [ITsi of effective area A ~ 10^ m^ is ~ 350 GeV. For the Lake Baikal detector fll 



the vertical muon threshold energy is ~ 150 GeV and the effective area is ~ 10^ m^. 
The signal in the detector, 5* = Nfj_{E^), must be compared with the square root of 
noise, a//V, given by 

y/N = ^Ax 6tx 1^(6) X 6e\ (9) 

where 5t is the observed duration of the burst, I^{0) is the background muon intensity 
(i.e., number of muons per unit area, time and solid angle) due to cosmic ray induced 



atmospheric showers [|1^, ^7j, and 69 is the angular resolution of the detector which 
is typically a few degrees. The smallest value oi n = S/ ^/N necessary to consider an 
excess signal as a positive detection depends on the search strategy of the particular 
detector. 
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TABLE I 

Expected number of muons in AMANDA detector for a GRB at 2; = 0.1 
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In Table I we present for illustration the result of our calculations of the number of 
muons expected in AMANDA and the corresponding signal to square root of noise 
ratios {S/\/N) due to a single GRB at redshift z = 0.1, for various possible values 
of its source luminosity (L), maximum energy cutoff at source (-E^max), and integral 
spectral index {a). The duration of the burst is taken to be 6t = Usee, somewhat 
intermediate between short- and long duration bursts. The angular resolution of the 
AMANDA detector has been taken to be 66 = 1°. The muon numbers are given for 
two different values of the zenith angle 6 of the source. The absorption of the TeV 
photons in the intergalactic infrared and optical backgrounds has been taken into 
account in calculating the numbers shown in Table I for which we have taken the 
optical depths corresponding to the higher level of the intergalactic infrared radiation 
field (IIRF) given in Ref. [0 in order to have conservative estimates of the number of 
muons produced. Throughout this paper we use optical depths corresponding to this 
higher level of IIRF whenever we include the effect of absorption of the TeV photons 
in the intergalactic space. 



Clearly, the number of muons ( "signal" ) and the signal-to-noise ratios {S/ y N) depend 
crucially on the various GRB parameters. Some of these dependencies are illustrated 
in more details in Figures below. As pointed out in |T^, for a given set of other 
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parameters, the smaller number of muons at the larger zenith angle is a consequence 
of the increase of the energy threshold with zenith angle because the muons have to 
traverse a greater amount of matter to reach the detector. Indeed, for too low value 
of the source gamma-ray energy cutoff -E-ymax, for sufficiently large zenith angle, the 
threshold gamma-ray energy E-yth required to produce muons above the threshold 
energy of the detector can be larger than -E^max, thus yielding no muons. Also, as 
expected, for a given luminosity L and spectral index a < 1, higher source cutoff 
energy -E^max yields higher number of muons as well as higher S/ \/N. However, the 
effect of decrement of the signal due to smaller -E^max can, in some cases, be over- 
compensated for by a harder source spectrum (i.e., smaller value of a) yielding a 
relatively larger signal. 

Fig. 1 and Fig. 2 show the number of muons expected in AMANDA as a function of 
gamma-ray integral spectral index of the GRB source assuming L = 3 x 10^^ ergs/ sec, 
Ste = 1 sec (so that Eqrb = 3 x 10^^ ergs), -E-ymax = lOTeV, and z = 0.1, for zenith 
angles 9 = 0° (Fig. 1) and 9 = 60° (Fig. 2). The results for the BAIKAL detector are 
shown in Fig. 3. Curves are shown for both without and with absorption of the TeV 



photons in the IIRF. In Fig. 1 and Fig. 2, the results of Ref. are also shown for 
comparison. 

From Fig. 1 and Fig. 2, we see that there are significant differences between our 



results and those of Ref. for the same set of GRB parameters, with our results 



for the number of muons being always consistently lower than those of Ref. ||T^. In 



the case of absence of absorption, we guess the difference between our results and 
those of [0 could be due to our using a specific cosmological model — as already 
mentioned, we use the "standard" cosmological model with a spatially fiat Universe 
with = 0.7 , ilm = 0.3 and Hq = 65 km sec~^ Mpc~^ — whereas Ref. does not 
explicitly mention the values of various cosmological parameters used. 

In the case when absorption of TeV photons in the IIRF is included, the difference 
between our results and those of is even more than that in the case of no absorp- 
tion. The flux of the TeV photons, and hence the number of muons in the detector, 
depend quite sensitively on the optical depth of the TeV photons, which in turn is 
governed by the strength of the IIRF. A relatively higher (lower) level of IIRF gives 
relatively larger (smaller) optical depth, resulting in relatively smaller (larger) num- 
ber of muons in the detector. Since the level of IIRF is not known with certainty, we 
have used optical depths corresponding to the higher level of IIRF from Ref. in 
order to obtain conservative estimates of the number of muons in the detector. How- 



ever, Ref. does not explicitly mention whether the optical depths corresponding 
to the higher or the lower level of the IIRF given in Ref. was used in their actual 
calculation, and so it is difficult to find the exact reasons for the differences between 
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our results and those of Ref. [0. For redshift z = 0.1 our results differ from the 
results of Ref. [0 by almost a factor of ten in the case when the effect of absorption 
of the TeV photons in the IIRF is included. 



The above discussion demonstrates the fact that the detectability of possible TeV 
gamma-rays from GRBs in ground based muon detectors depends sensitively on the 
various cosmological parameters and the strength and spectrum of the IIRF in addi- 
tion to the intrinsic GRB parameters. Conversely, within the context of the "standard 
cosmological model" , the detection or non-detection of TeV photons from GRBs can 
significantly constrain the strength and spectrum of the IIRF, provided the intrinsic 
GRB parameters such as their total luminosity, redshift and duration are known with 
reasonable accuracy from independent observations. 

In Fig. 4. the number of muons expected in AMANDA from a GRB is shown as a 
function of the redshift of the GRB for -E^max = 10 TeV and for two different cases, 
namely, L = 3 x 10^^ ergs/ sec , a = 1.0, and L = 10^^ ergs/ sec, a = 0.5. The 
absorption of TeV photons is included with optical depth corresponding to the higher 
level of IIRF from 0. As expected, the number of muons from a GRB falls steeply 
with increasing redshift due to absorption of TeV photons, and so detection is possible 
only for relatively nearby (low redshift) bursts. 

Fig. 5. shows the expected number of muons (for 6 = 0°) in AMANDA as a function 
of the total luminosity of the GRB for a fixed duration St = Isec of the burst. In 
Figs. 6-8 we show the effect of varying the physical parameters of a GRB on the signal 
to square root of noise ratios in the AMANDA detector. It is clear that high signal 
to square root of noise detection of TeV photons from GRBs in AMANDA (or other 
ground based muon detectors of similar area and energy threshold) will be possible 
mostly for relatively nearby {z < 0.1), long duration ( > several tens of seconds), 
high luminosity ( > few x 10^^ ergs/ sec) and hard spectrum [a < 1) GRBs. 



Table II 



BATSE 
Number 


50-300 keV Luminosity 
(Ll) from |20| (ergs/sec) 


Redshift 

z 


Burst duration 
{6t) (sec) 


1 MeV-10 TeV Luminosity 
(Lh) (ergs/sec) 


6707 


2.63 X 10^^ 


0.0085 


20.6 


2.48 X 10^^ 


2123 


0.18 X 10^° 


0.1 


22.0 


10.06 X 10^^ 


2316 


0.64 X 10^^ 


0.1 


29.2 


7.58 X 10^^ 


3055 


0.61 X 10^° 


0.2 


40.6 


2.12 X 10^^ 


3352 


0.62 X 10^° 


0.2 


46.3 


1.86 X 10^^ 


4368 


6.04 X 10^^ 


0.4 


36.5 


0.88 X 10^° 



10 



In Table II we list some of the nearby GRBs from the BATSE catalogue whose lumi- 
nosities in the 50 to 300 keV photon energy band are known and whose redshifts have 
been estimated in Ref. by using a conjectured variability - luminosity relationship 
for GRBs which is calibrated using seven GRBs whose redshifts are known indepen- 
dently from afterglow studies. We have calculated the luminosities for photon energies 
in the range 1 MeV to 10 TeV required for these GRBs to produce S/\/N = 1.5 at 
zenith angle 6* = 0° in AMANDA. We assume a = 1.0 for these GRBs. The luminosi- 
ties so calculated, denoted by L^, are displayed in the last column of Table II. The 
intergalactic absorption of TeV photons corresponding to the higher IIRF level from 
has been taken into account in this calculation. 

From Table II we see that the high energy (1 MeV to 10 TeV) luminosities (Lh) 
required to detect TeV photons from these GRBs in AMANDA are much higher than 
their luminosities Ll in the lower energy (50 keV to 300 keV) band. No TeV photons 
from these GRBs have been reported to be detected by AMANDA. The most likely 
explanation is that these GRBs did not emit sufficient number of TeV photons that 
would create sufficient number of muons above the threshold of AMANDA and/or 
the intergalatic absorption of TeV photons is even stronger than what we assumed. 

4 Average Rate of Detectable TeV GRBs in Muon 
Detectors 

In the previous sections we have studied the detectability of TeV photons from indi- 
vidual GRBs in muon detectors by calculating the number of muons produced and 
the signal to noise ratios as functions of various intrinsic parameters of the bursts 
such as their luminosities, integral spectral index, redshift, and burst duration. 

In this section we estimate the expected rate of positive detection of TeV photons 
from all GRBs in the Universe in muon detectors such as AMANDA. This requires the 
knowledge of the forms of the luminosity function (LP) and the redshift distribution 
(RD) of the GRBs in the Universe. Unfortunately, direct measurements of redshifts 
exist for only a few GRBs (from optical afterglow measurements), and so the RD and 
the LF of the GRB population are not known with certainty. Nevertheless, recently 
there have been attempts to indirectly derive the redshifts (and hence luminosity) 
distributions of GRBs from the BATSE data exploiting various features of the light 
curves of individual GRBs in the BATSE data. 
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4.1 Luminosity Function and Redshift Distribution of GRBs 



Recently, authors of Ref. have found that the luminosities of seven GRBs in 
the BATSE catalogue with known redshifts (measured from their optical afterglows) 
are correlated (as a power law) with the variability of the individual bursts, the 
variability being defined as the normalized variance of the observed 50-300 keV light 
curve about a fitted smooth light curve. Based on this observation, the authors 
of Ref. [^] have conjectured a variability (\^)-luminosity (L) relationship which is 
hypothesized to be true for all GRBs in the BATSE catalogue. Use of this conjectured 
V-L relationship, which is calibrated using the seven GRBs with known redshifts, 
then allowed estimation of the redshifts of about 220 GRBs in the BATSE catalogue, 
which in turn allowed construction of a LP and a RD for GRBs in general. 



Independently, authors of Ref. |^ have found a correlation between luminosity and 
the so-called "lag" (xiag), which is the time delay between the peaks in the light 
curves of individual GRBs. It was then pointed out in Ref. that the above two 
relationships, namely the V-L and the riag-L relationships, if true in general, would 
together imply a V-Ti^g relationship which could be tested directly on the measured 
data for BATSE observed GRBs without reference to their redshifts. Indeed, Ref. p2| 



claims strong concordance of the above mentioned ^-riag relationship with the ob- 
served BATSE data for 112 GRBs, based on which a LP has been derived for GRBs 
in general, which is a broken power-law with a break at L ^ 2 x 10^^ erg/ sec, with 
LP going as /.-I 'i'iO-i below the break and /,-2.8±o.2 qJ^q^q break. The RD giving 



the number density of GRBs as a function of redshift is also derived [^], which goes 

as (1 + 2)2-5±0.3_ 

We must remember that the LP mentioned above refers to the luminosity in the 
BATSE band (50-300 keV), which we have denoted above by Ll. We are, however, 
interested in the LP expressed as a function of the luminosity in the 1 MeV — 10 TeV 
interval, which we have denoted above by Lh- It is, however, easy to see that, for a 
single power-law photon spectrum across all energies, which we assume for simplicity 
in this paper, the LP expressed as a function of Lh has the same functional form as 
that in terms of Ll, except that the break of the LP occurs at a different value of 
luminosity which depends on the integral spectral index a. For a given GRB with a 
definite Ll and Lh, one can see that the two are related as 

Lh _ (10TeV)-"+^-(lMeV)---+^ _ 

Ll (300keV)-"+i-(50keV)-°+i ^ ^' ^ ' 

which is valid for all GRBs assumed to have the same integral spectral index a. 
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Based on the above considerations, we shall assume the following LF for GRBs: 




for Lh < Lh* 
for Lff > Lh* 



where the break luminosity Lh* is given by Lh* = C{a)LL*, with Ll* 
10^^ erg/ sec being the break of the LF derived in Ref. p2| . 



fill 



2 X 



Below, we shall often drop the subscript H in Lh, and so, unless otherwise specified, 
L will mean Lh- 

For the rate-density of GRBs (number of bursts per unit volume and unit time) in the 
Universe as a function of redshift, we will assume it to be proportional to (1 + z)"^'^ 
for the redshift range of our interest (up to z ~ 4). This redshift distribution 
follows star formation rate (SFR) ||23| for z < 2, but unlike the observed SFR based 
on optical observations, which generally yield an SFR that either levels off or falls 
with redshift for z > 2, the above ORB rate increase monotonically with z at least 
up to z ~ 5 



2^ and possibly even up to 2; ~ 10 pOf- The suggestion pO|, |22 



IS 

that the SFR based on optical observations may have underestimated the true SFR 
at high redshifts because of cosmological reddening effects, and since gamma rays 
are not subject to the reddening effect, the GRB rate may indeed reflect the true 
SFR at high redshifts provided, of course, the physical connection between the GRB 
phenomenon and the star formation process is clearly understood by future studies. 

With the LF and rate-density (i.e., redshift distribution) of GRBs specified above, we 
now proceed to calculate the expected average rate of detection of TeV GRBs in the 
muon detectors. In doing this we follow the formalism described in Ref. p4 . 



4.2 Calculation of the Rate of TeV GRBs 

The first step is to calculate the redshift up to which a GRB of a given L, E^-^ax, « 
and St will be detectable at a zenith angle ^ in a given detector. By detectable we 
mean the statistical significance of the signal, i.e., the signal to square root of noise 
ratio, S/ \/N, is larger than some preassigned value. 
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TABLE III 



Luminosity (L) 
(ergs/sec) 


Redshift up to which 
the GRB is detectable 


10^1 


0.057 


1052 


0.091 


1053 


0.129 


1054 


0.168 


1055 


0.207 


1056 


0.245 



Table III illustrates the values of maximum redshifts up to which GRBs of various 
L are detectable by AMANDA at zenith angle ^ = 0° with S/\fN > 1.5, assuming 
E-yrnax = 10 TcV, « = 0.8 and 6t = 20 sec, with absorption of TeV photons taken into 
account. Fig. 9 shows this in more details for three different values of the integral 
spectral index a. The regions below the curves in this figure are the allowed regions 
for positive detections with S/\fN > 1.5. 

Next, we calculate the fraction, f{L,cos6), of all bursts of given luminosity L that 
are detectable at a zenith angle 6 in the sense defined above. To do this, we first 
define, for a burst of luminosity L at redshift z at a zenith angle 6 with respect to 
a given detector, a quantity J(L, z, cos 9) such that J(L, z, cos ^) = 1 if the burst is 
detectable and J{L, z,cos6) = if it is not. With the redshift distribution of GRBs 
as specified above, we can write 

f(r - Io-JiL,z,cose)il + zr%dz 

where Zmax is the maximum redshift of the GRBs in the Universe, and dV{z) is the 
volume element of the Universe between redshifts z and z + dz, which is given by 



dV 47r (i r , , 1 , , 

dz [l + zfdz^ ^ " ^ ^ ' 

In the above expression r{z) is the comoving distance as defined in Eq. (^. 

In our numerical calculations we shall take Zmax = 4. 



Finally, folding the fraction / with the LF given in Eq. (|TT]) and the zenith angle 
distribution of the bursts, dS {cos 9) /d cos 9, we write the expected TeV GRB rate in 
muon detector, -Rxev, as 
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i?TeV = i?total / / — rf^/ COS ^ — ^ -^dLd(cOS 6) , 14 

JLmin io.5 aL a COS 6' 

where dN{L)/dL is now the normahzed LF and i?totai is the total observed GRB rate, 
and we assume that the detector can view GRBs with good efficiency between zenith 
angles ^ = 0° and 9 — 60°. We shall assume a uniform zenith angle distribution for 
the bursts and set dS {cos 9) /d cos 9 — 1/2. Also we take Lmin = 10^^ erg/ sec and 
i-max = 10^^ erg/ sec. 

BATSE observes about 1 burst per day with a detection efficiency of 0.3. Thus, the 

GRB rate deduced from BATSE observations is Rb ~ 1000 yr~^. Assuming that 
BATSE can observe GRBs up to a maximum redshift Zmax = 4 we set i?totai = -Rb- 

We have calculated -RtcV for AMANDA for a wide range of values of the intrinsic 
parameters specifying the bursts. These rates turn out to be rather low. For example, 
with £^-ymax = lOTeV, a = 0.8 and dt = 20 sec, and requiring S/^/N > 1.5 (the same 
parameter set as used in obtaining the numbers in Table III), we get ~ 1 detectable 
TeV GRB in 25 years in AMANDA. The rate depends somewhat on a; for example, 
for a = 0.5 with other parameters kept same as above, the rate improves to ~ 1 burst 
in 20 years, while for a = 1, the rate is ~ 1 in 40 years. In ICECUBE, assuming 
its effective area will be about 10 times larger, the above expected rates will be 
correspondingly larger by the same factor, assuming same threshold energy for muon 
detection as in AMANDA. 

Note also that the above rates were calculated for only a 1.5cr detectability; rates for 
higher-(7 detectability will be even lower. 

4.3 Discussion 

The reasons for the rather low detectability rate of TeV photons from GRBs in existing 
muon detectors are clear: The absorption of TeV photons in the inter galactic infrared 
background allows detection of TeV photons from only relatively nearby {z < 0.1) 
bursts as clear from Fig. 9, unless the luminosity of the burst is relatively large 
( ^ 10^^ erg/ sec). However, for a LF falling with L, and a redshift distribution scaling 
like the SFR (which increases with z), low-redshift and high- luminosity GRBs arc 
relatively rare. Nevertheless, as we have demonstrated in Figs. 6-9, occasional nearby 
GRBs of sufficiently long duration, large luminosity and hard spectrum are likely to 
be detectable in existing detectors like AMANDA, and will certainly be detectable in 
the next generation detectors such as ICECUBE. 
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5 Summary and Conclusions 



It is possible that GRBs emit not only sub-MeV photons as detected in satellite- 
borne detectors, but also higher energy photons extending to TeV energies. For a 
GRB photon spectrum falling with energy, as is usually the case, the non-detection 
of TeV photons in satellite-borne detectors could be due to their limited size. On 
the other hand, for sufficiently hard photon spectrum with integral spectral index 
a < 1, although the total number of photons is dominated by those at the low (keV 
- MeV) energy end, the total energy of the burst would be carried away mostly by 
the few high (say, TeV)- energy photons. Thus, the actual total energy emitted in a 
burst may be much larger than the total energy estimated from the observed fluence 
of sub-MeV photons in the satellite detectors, which, if true, would have tremendous 
implications for the theories of origin of GRBs. It is, therefore, very important to 
study possible ways of detecting possible TeV photons from GRBs. 

Refs. [0, |17| suggested that TeV photons from GRBs might be detectable in the 
existing muon detectors such as AMANDA, Lake Baikal, and Milagro and future de- 
tectors such as the proposed ICECUBE, by detecting the muons in the atmospheric 



showers created by TeV photons. Following the calculations of Ref. IT^, we have 
made detailed calculations of the detectability of possible TeV photons from individ- 
ual GRBs by AMANDA-type detectors, as functions of various GRB parameters such 
as their luminosity, duration, integral spectral index, and redshift. While our calcula- 
tions generally yield lower number of expected muons per burst than that calculated 
m Ref. in for the same set of GRB parameters, we conclude that sufficiently high 



luminosity, long duration, and hard spectrum nearby GRBs would still be detectable 
in AMANDA-class detectors and will certainly be detectable in next generation de- 
tectors such as ICECUBE. We have also estimated the expected rate of TeV GRB 
events in these detectors using recent information on the luminosity function of GRBs 
and their redshift distribution in the Universe. Our conservative estimates, including 
the effect of absorption of TeV photons in the intergalactic space, show that while 
the expected average rate of TeV GRB events in AMANDA is rather low — about 
one GRB in 20 years for sufficiently hard (integral spectral index a = 0.5) spectrum, 
the rate in up-coming bigger detectors such as ICECUBE could be about 10 times 
larger, i.e., about 1 in every couple of years. 
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Figure 1: Number of muons per burst at zenith angle ^ = 0° in AMANDA as a 
function of the integral spectral index of the GRB photon spectrum, for a GRB at 
redshift z = 0.1, total energy Eqjiq = 3 x 10^^ erg, and maximum cutoff energy 
E^raax = 10 TcV. Curvcs are shown for with and without absorption of the photons 
in the intergalactic infrared field (IIRF). Corresponding curves from the calculation 
of Ref. |T^, indicated above by A H (their Fig. 1) for the same set of parameter values 
as above, are also shown for comparison. In the case of absorption, our results are 
almost a factor of 10 below those in A H (see text for discussion). 
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Figure 2: Same as Fig. 1, but for zenith angle 9 = 60° 
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Figure 3: Number of muons per burst in the Lake Baikal detector as a function of 
the integral spectral index of the GRB photon spectrum, for the same set of GRB 
parameters as in Fig. 1. Curves are shown for with and without absorption of the 
photons in the IIRF and for two different zenith angles ^ = 0° and 9 — 60°. 
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Figure 4: Number of muons per burst in the AMANDA detector at zenith angle 
^ = 0° as a function of the redshift of the burst for various GRB parameters as 
indicated, taking into account the absorption of photons in the IIRF corresponding 
to the higher level of IIRF (see text). 
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Figure 5: Number of muons per burst in the AMANDA detector at zenith angle ^ = 0° 
as a function of the luminosity of the burst for various GRB parameters as indicated. 
Curves are shown for both with absorption (dashed curves) and without absorption 
(solid curves) in the IIRF. The upper solid or dashed curves are for maximum spectral 



cutoff of the source at E, 



7max 



10 TeV and the lower ones for E. 



-ymax 



5 TeV. 
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Figure 6: The signal to square root of noise ratio in the AMANDA detector as a 
function of the integral spectral index of the GRB photon spectrum, corresponding 
to the muon numbers shown in Figs. 1 and 2. 
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Figure 7: The signal to square root of noise ratio in the AMANDA detector as a 
function of redshift of the GRB, corresponding to the muon numbers shown in Fig. 4. 
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Figure 8: The signal to square root of noise ratio in the AMANDA detector as a 
function of the duration of the GRB. 
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Figure 9: The maximum redshift of a GRB that can be detected with signal to square 
root of noise ratio S/^/N > 1.5 in AMANDA as a function of the luminosity of the 
burst for several choices of integral spectral index a and for the indicated values of 
-E-ymax and burst duration. The regions below the respective curves are the allowed 
regions for positive detection. 
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